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ABSTRACT 

N2H + observations of molecular cloud cores in Taurus with the Nobeyama 45 
m radio telescope are reported. We compare "cores with young stars" with "cores 
without young stars". The differences in core radius, linewidth, and core mass 
are small. Linewidth is dominated by thermal motions in both cases. N2H + maps 
show that the intensity distribution does not differ much between cores without 
stars and those with stars. This is in contrast to the result previously obtained 
in H 13 CO + toward Taurus molecular cloud cores. Larger degree of depletion of 
H 13 CO + in starless cores will be one possible explanation for this difference. We 
studied the physical state of molecular cloud cores in terms of "critical pressure" 
for the surface (external) pressure. There is no systematic difference between 
starless cores and cores with stars in this analysis. Both are not far from the 
critical state for pressure equilibrium. We suggest that molecular cloud cores 
in which thermal support is dominated evolve toward star formation by keeping 
close to the critical state. This result is in contrast with that obtained in the 
intermediate-mass star forming region OMC-2/3, where molecular cloud cores 
evolve by decreasing the critical pressure appreciably. We investigate the radial 
distribution of the integrated intensity. Cores with stars are found to have shallow 
(—1.8 to —1.6) power-law density profiles. 

Subject headings: ISM: clouds — ISM: individual (L1495, L1521F, L1527, L1551, 
Taurus Cloud Complex) — ISM: molecules — ISM: structure — radio lines: ISM— 
stars: formation 
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1. Introduction 

It is of great interest when and how the process of star formation takes place. Nakano 
(1998) theoretically studied the condition for the onset of star formation on the basis of the 
virial theorem. For the surface (external) pressure P s , we can define the critical pressure 
P cr , which is a function of the effective sound speed and the mass of the molecular cloud 
core (Spitzer 1968; Nakano 1998). There are two equilibria if P cr > P s , whereas there is 
no equilibrium if P cr < P s . Nakano (1998) suggested that dissipation of turbulence in the 
core will result in the decrease in P cr , which will lead to the onset of star formation because 
of the absence of equilibrium state. Aso et al. (2000) observed the intermediate-mass star 
forming region OMC-2/3 in Orion by using the H 13 CO + J = 1 — > line. They found that 
cores without young stars tend to have P cr > P s and that cores with young stars tend to 
have P cr < P s . This fact suggests that dissipation of turbulence leads to the onset of star 
formation in OMC-2/3. 

It seems that molecular cloud cores in low-mass star forming regions like Taurus have 
lower surface (external) pressure than those in massive or intermediate-mass star forming 
regions like Orion(Tatematsu et al. 1993). The present study is aimed to investigate whether 
dissipation of turbulence plays an important role also in low-mass star forming regions. In 
low-mass star forming regions, the choice of the molecular-line tracers is crucial. It is known 
that molecular depletion seems to be serious for CO, HCO + , CCS, CS, and their isotopomers 
in starless, cold molecular cloud cores (Caselli et al. 1999; Aikawa et al. 2001; Bergin et al. 
2001, 2002; Li et al. 2002; Tafalla et al. 2002; Lee et al. 2003; Shematovich et al. 2003). On the 
other hand, N-bearing molecules, NH 3 and N 2 H + do not show strong evidence of depletion, 
except for the very center of the cores (Bergin et al. 2001, 2002; Caselli et al. 2002b; Tafalla 
et al. 2002; Lee et al. 2003; Shinnaga et al. 2003). Furthermore, although H 13 CO + is one of 
the popular tracers of the molecular cloud core, it is known that the linewidth of the H 13 CO + 
core can be broaden by the molecular outflow, and this is more serious in studying low-mass 
star forming regions (Aso et al. 2000). We employ N 2 H + , which is less affected by depletion 
or molecular outflow. N 2 H + is found to trace the quiescent molecular gas (Womack et al. 
1993; Bachiller & Perez Gutierrez 1997), and this choice would minimize the influence of the 
molecular outflow. 

The target molecular cloud cores are those observed by Mizuno et al. (1994) in H 13 CO + . 
Molecular cloud cores associated with young stars are called "cores with stars." Molecular 
cloud cores not associated with young stars are called "starless cores." Our samples contain 
three starless cores and four cores with stars. We happen to detect one new starless core 
during the observations. The distances to cores in Taurus are assumed to be 140 pc. At this 
distance, 1 arcmin corresponds to 0.041 pc. 
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2. Observations 

We carried out our observations by using the 45 m radio telescope of Nobeyama Radio 
Observatory from 2003 January 7 to 9. The receiver front end was the 25-element double- 
sideband focal-plane SIS receiver BEARS, which has a beam separation of 41'/1. We observed 
molecular cloud cores in N 2 H+ J = 1 -> at 93.1737767 GHz (Caselli et al. 1995). We 
mapped regions by shifting the telescope by half the beam separation, so the final map 
consists of data on a 20'.'55 spacing grid. During the observations, one element of BEARS 
was not available due to a technical problem, so we used 24 elements. The half-power 
beamwidth of the element beams of BEARS was estimated to be 17 / /8±0 / . / 4 at 93 GHz, 
which was close to the data grid spacing. The employed receiver back end was a digital 
autocorrelator. We selected the spectral resolution to be 37.8 kHz (corresponding to ~ 0.12 
km s -1 ). Spectra were obtained in the position-switching mode. To compensate for the daily 
intensity scale variation, we observed a local intensity peak at R.A. (J2000) = 5 ft 35 m 16!0, 
Decl.(J2000) = -5°19'26" near Orion KL with BEARS every day. Furthermore, we observed 
this intensity peak using a single-sideband SIS receiver S100 to calibrate the gain of each 
BEARS element and to establish the absolute intensity scale. The line intensity is reported 
in terms of the corrected antenna temperature T\ to be observed with the S100 receiver 
in this paper. The main-beam efficiency with S100 at 93 GHz is estimated to be 0.515, 
by interpolating the measurements at 86 and 100 GHz by the observatory. The telescope 
pointing was established by observing the SiO maser source NML Tau at 43 GHz every 
1 — 1.5 hours during the observations. The data were reduced by using the software package 
NewStar of Nobeyama Radio Observatory and IDL of Research Systems, Inc. 

Prior to our observations, sixty-three low-mass molecular cloud cores have been observed 
in N 2 H + by Caselli et al. (2002c), but with lower angular resolution (54 arcsec). We will 
refer to their results to secure our discussion on the basis of our observations toward the 
eight cores but with better angular resolution (18" 6 beam and 20'.'55 grid). Since the typical 
radius of the molecular cloud core is about 50 arcsec in the region, the improvement in spatial 
resolution helps us to derive the physical parameters more precisely. 

3. Results 

Figure 1 shows that the distribution of the velocity-integrated intensity of the main 
N 2 H + J = 1 — > component (Fi, F = 2, 3— >1, 2) toward the eight molecular cloud cores. 
The core associated with L1551 IRS5 (L1551) is 2—3 times as intense as the other cores, and 
we double the contour interval for this source for clarity. We use the main component rather 
than the optically thinner components to obtain better signal-to-noise ratios. For checking, 
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we have also made the integrated intensity map of the Fi, F — 0, 1-^1, 2 component (not 
shown, the map quality is much worse), and have confirmed that the intensity peak position 
is consistent in general. As shown later, the optical depth is not very large even for the main 
component (~ unity). The positions of the intensity maxima of the cores are summarized 
in Table 1. The core name with the prefix of Miz refers to the core number in Mizuno et al. 
(1994). Core Miz8b is newly found in our observations. Core L1527 is associated with a Class 
protostar, while cores Miz7, Miz8, and, L1551 are associated with Class I protostars. These 
are "cores with stars." Furthermore, a Class protostar L1551 NE (Barsony & Chandler 
1993) is located near the core boundary of L1551, but its relationship to the N 2 H + core is 
not clear in our map (see Saito et al. (2001); Yokogawa et al. (2003) for a detailed study of 
this region in C 18 and CS). Cores Mizl, Miz2, L1521F, and Miz8b are "starless cores". 

The basic physical parameters of the cores are summarized in Table 2. The HWHM core 
radius R is measured as y/S/ ir (S is the core area S at the half maximum), and then corrected 
for the telescope beam size. There is no remarkable difference in core radius between starless 
cores (R = 0.035±0.004 pc) and core with stars (R = 0.031±0.006 pc). The study of 
Caselli et al. (2002c) shows that cores with stars are slightly larger than starless cores, but 
difference is not large. Our result is in marked contrast with that of Mizuno et al. (1994) 
in H 13 CO + , in which cores with stars tend to be more compact. Depletion of H 13 CO + will 
be a plausible explanation, because it is known that N 2 H + is more robust for depletion than 
H 13 CO+ (Bergin et al. 2001; Caselli et al. 2002b; Lee et al. 2003). It is suggested that 
depletion becomes less prominent in cores with stars because molecules evaporate from the 
dust grain due to radiation and outflow from the newly formed star (Langer et al. 2000). 
If some molecule is depleted in starless cores, we expect that intensity tends to be weak 
in starless cores, and the radial intensity profile is more flat-topped (because depletion is 
more effective at the high-density core center), causing larger radius in starless cores. The 
tendency observed in H 13 CO + can be explained if this molecule is depleted substantially only 
in starless cores. 

Next, we derive the line optical depth, linewidth, and mass. We fit the observed spec- 
trum by using the hyperfine spectrum model consisting of multiple Gaussian components 
with line optical depth effect. The intrinsic relative intensities of the hyperfine components 
are taken from Tine et al. (2000). The free parameters are the excitation temperature T eX) 
the sum of optical depths of the hyperfine components ttot, systemic velocity (radial veloc- 
ity), and intrinsic linewidth (which is not broadened due to line optical depth). The details 
of the column density estimation are given in Caselli et al. (2002b). We used the optically 
thick formula for the central 3x3 positions, and the optically thin formula for the weaker, 
outer part. Figure 2 shows the results of the hyperfine fitting, and Table 2 lists the physical 
parameters from the hyperfine fitting. The excitation temperature at the intensity peak is 
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T ex = 6.2 ± 2.3 K for the seven cores excluding Mizl. That for the composite spectrum from 
the central nine (3x3) positions is T ex = 5.7 ±1.2 K for the eight cores. The optical depth 
of the main component (F 1 , F = 2, 3— »1, 2), which is equal to 0.259 ttot, is found to be 
moderate (1.54±0.78 for the intensity peak and 1.07±0.65 for the central nine positions). 
The derived linewidth is then corrected for the frequency resolution of the spectrometer. 
Figure 3 plots the intrinsic linewidth against the angular distance from the core center (im- 
pact parameter) b. This corresponds to the Type 4 (single-tracer, single-cloud) linewidth-size 
relation in Goodman et al. (1998). The power- law index of the linewidth-size relation ob- 
tained by non-linear least-squares fitting is listed in Table 3. Caselli et al. (2002c) showed 
that molecular cloud cores have a variety of linewidth-size relation; cores show positive, flat, 
and negative correlation of linewidth with the impact parameter. Our samples show that 
linewidth decreases or constant with increasing b in general (Table 2). The larger linewidth 
at the center in cores with stars could be due to the influence of protostellar collapse (Zhou 
et al. 1994; Caselli et al. 2002a) and/or molecular outflow (Aso et al. 2000), although we 
selected the molecular line to minimize such effects. To discuss the intrinsic core properties 
by eliminating this effect, we use the average linewidth for b > 40" in Miz8 and L1527 and 
for b > 80" in L1551. The intrinsic N 2 H+ linewidth of the starless core is Av = 0.256±0.024 
km s -1 , while that of the core with star is Av = 0.309±0.070 km s _1 . So, there is no signif- 
icant difference. By assuming a gas kinetic temperature of 10 K, we derive the non-thermal 
linewidth Avnt and the total linewidth Avtot (Fuller & Myers 1992). Avnt is estimated 
to be 0.224±0.030 km s _1 and 0.284±0.082 km s -1 for starless cores and cores with stars, 
respectively. Avtot is derived to be 0.497±0.013 km s -1 and 0.528±0.042 km s" 1 for starless 
cores and cores with stars, respectively. The thermal linewidth Avt for the mean molecular 
weight (2.33 m H ) at 10 K is 0.443 km s _1 . So, both starless cores and cores with stars are 
dominated by the thermal support. The effective pressure is 1.26±0.07 and 1.43±0.23 times 
as large as the thermal pressure for starless cores and cores with stars, respectively. Table 3 
lists the H 2 column density iV(H 2 ) obtained for the spectrum at the intensity peak, and for 
the composite spectrum from the nine (3x3) positions centered on the intensity peak. We 
derive the core mass M by integrating the column density over the core. We adopted the 
abundance of N 2 H+ to be 3.0xl0" 10 relative to H 2 (Caselli et al. 2002c). The core mass M 
is listed in Table 3. The average value is M = 1.30±0.63 M & for the starless cores, and is M 
= 1.59±0.97 M for the cores with stars. These values are not much different, taking into 
account the uncertainty in mass estimation (a factor of 2). 
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4. Discussion 

We investigate whether the cores have equilibrium states or not by following the discus- 
sion of Nakano (1998) (see §1 of the present paper for a brief summary). The virial theorem 
is useful not only for stable states but also for the unstable states not very far from the 
critical state (Nakano 1998). For consistency, the formula in Aso et al. (2000) is used again, 
P cr = l/12nG 3 M 2 x (5/3) 3 x (9/4C e 2 // ) 4 ! where the effective sound speed C eff is C eff = 
(Av% OT /8 In 2) 1 / 2 . Table 3 lists the critical pressure P cr . P cr is derived to be (5.4±2.6)xl0 5 
and (5.9±2.2)xl0 5 K cm -3 for starless cores and cores with stars, respectively, when we 
do not correct for the decrease in mass (see below). When we correct for the decrease in 
mass due to protostellar collapse and outflow for cores with stars, the critical pressure P cr 
is derived to be (5.4±2.6)xl0 5 and (3.5±1.6)xl0 5 K cm -3 for starless cores and cores with 
stars, respectively. We estimate P s in Taurus to be of order 4xl0 5 K cm -3 from the col- 
umn density of low-density molecular gas and the coefficient of the linewidth-size relation 
(Tatematsu et al. 1993; Onishi et al. 2002). This means that both starless cores and cores 
with stars are close to the critical state P cr ~ P s in Taurus. 

To know whether turbulence dissipation leads to the onset of star formation, we need 
to know the physical condition prior to star formation. The derived mass of the molecular 
cloud core is the current value. This may have decreased from the original value, because 
part of the original core mass has lost through star formation. We estimate the decrease in 
the core mass due to these star formation activities as follows. Part of the core mass will 
be accreted onto the protostar, and part of it will be swept up by the molecular outflow. 
The high-velocity wind from the protostar will sweep up the ambient matter (mostly the 
parent molecular cloud core, and possibly also the less dense envelope), and this entrained 
mass together with the wind mass constitutes to the outflow mass. The mass of the T Tauri 
star is typically 0.5 M Q . The Class I protostar is supposed to have accreted more than half 
of the final stellar mass (Bontemps et al. 1996). The mass of the forming star is estimated 
to be 0.1-0.15 M for IRAS 04169+2702 from the accretion luminosity estimate (Ohashi 
et al. 1997b). The mass of the central star of L1551 IRS5 is estimated to be 0.15—0.5 
M & from the kinematics of the protoplanetary disk (Momose et al. 1998) and to be 0.7—3 
M & from the stellar luminosity and disk mass (Saito et al. 1996). The mass of the central 
star of IRAS 04368+2557 (L1527) is estimated to be 0.1 M Q from the kinematics of the 
protoplanetary disk and the accretion luminosity (Ohashi et al. 1997a). Hogerheijde et al. 
(1998) estimated the maximum central star mass to be 2.6 and 0.2 M & for L1551 IRS5 and 
IRAS 04368+2557 (L1527), respectively, from the accretion luminosity. IRAS 04166+2706 
(Miz7), IRAS 04169+2702 (Miz8) and IRAS 04368+2557 (L1527) are known to accompany 
molecular outflows (Bontemps et al. 1996). The mass of molecular outflow associated with 
L1551 IRS5 is 0.3 M (Lada 1985). Hogerheijde et al. (1998) estimated the outflow mass 



-7- 



from 12 CO (3->2) to be 3.1 and 0.18 M for L1551 IRS5 and IRAS 04368+2557 (L1527), 
respectively. Taking them into account, we estimate the mass lost from the parent molecular 
cloud core is about 0.4, 0.4, 1.0, and 0.2 M for IRAS 04166+2706 (Miz7), IRAS 04169+2702 
(Miz8), L1551 IRS5, and IRAS 04368+2557 (L1527), respectively. 

Figure 4 plots the critical pressure P cr against the core mass M. There is no systematic 
difference between starless cores and cores with stars in Taurus. Regardless of starless cores 
or cores with stars, six out of the eight molecular cloud cores are located near the critical 
state of equilibrium. Two cores, L1521F and L1551 are located slightly below P s . However, 
P s may have local variation and mass estimation will be uncertain by a factor of 2. It is hard 
to conclude that these two cores are far from the critical state. Although there is a possibility 
that the core in Taurus slightly dissipates turbulence resulting in star formation and then 
the star formation activity increases the nonthermal linewidth, the present observations do 
not provide us with the evidence for the dissipation of turbulence. Within the accuracy of 
the current study, we conclude that both starless cores and cores with stars are close to 
the critical state. Myers & Benson (1983) obtained similar results from NH3 observations 
of low-mass molecular cloud cores; cores are close to the critical state for equilibrium and 
stability if the Doppler linewidth supports cores. Caselli et al. (2002c) studied the ratio of 
the core mass to the virial mass, and derived to be 1.3 ± 0.3 and 1.4 ± 0.3 for starless cores 
and cores with stars, respectively. These values are almost identical: both cores are close to 
virial equilibrium. We suggest that the thermally-supported cores evolve keeping the critical 
state for equilibrium. This result for the low-mass star forming region is in contrast with 
that obtained in the intermediate-mass star forming region OMC-2/3. In OMC-2/3, starless 
cores and cores with stars show clearly different states. 

We should revisit OMC-2/3 to see whether the results reported are still correct even 
if we take into account the depletion of H 13 CO + . The molecular cloud cores in OMC- 
2/3 is warmer (~ 20 K, Cesaroni & Wilson (1994)), and show active star formation. The 
sublimation temperature of the CO ice is 16 K in cloud cores (Langer et al. 2000). So, 
depletion will not be serious even for H 13 CO + . Aso et al. (2000) have not taken into account 
the mass accreted onto the protostar or swept by the outflow, but this correction will simply 
enhance the observed difference between starless cores and cores with stars. We conclude 
that the results obtained by Aso et al. (2000) is unchanged even in the present context. 

Lastly, we investigate the radial density profile on the basis of the radial distribution 
of the integrated intensity of the main component (Fig. 5). We use here the integrated 
intensity, because it is straightforward and more reliable than the hyperfine fitted column 
density for outer regions. The intensity profile is fitted with a power-law / oc r p convolved 
with the telescope beam. The fitting results are listed in Table 3. Mizl and Miz2 may have 
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(a hint of) central flattening. For Mizl and Miz2, the fitting result only for the outer part 
(> 35 arcsec or > 4800 AU) is shown as broken lines in Figure 5. The central flattening 
in cores in low-mass star forming regions was previously reported by Ward-Thompson et 
al. (1999); Andre et al. (1993); Caselli et al. (2002c) from the dust continuum observations 
and N2H + observations. However, the integrated intensity profile in Mizl and Miz2 can 
also be fitted reasonably with single (shallow) power laws (solid lines). The data quality of 
the present observation is not enough to distinguish these models. The power-law index of 
the radial distribution is p = —1.00 ± 0.14 (using the outer part for Mizl and Miz2) and 
—0.72 ± 0.07 for starless cores and cores with stars, respectively. When we assume that 
the column density is proportional to the integrated intensity, the power p in the intensity 
profile is related to the power a in the density profile p oc r a as « = p - 1. L1521F was 
observed in N 2 H + with the Berkeley-Illinois-Maryland Association (BIMA) millimeter array 
and with SCUBA on the James Clerk Maxwell Telescope (JCMT) (Shinnaga et al. 2003). 
These observations show the centrally peaked intensity profile of L1521F (from their Figure 
3, the BIMA intensity profile is fitted with p ~ —1.3 and the SCUBA intensity profile is 
fitted with p ~ — 1.0 for r > 10 arcsec). These data do not show the central flattening for 
L1521F, when we take into account their map resolution. Some starless cores (Mizl and 
Miz2) may have central flattening, while some starless cores may not. On the other hand, 
there is no hint of central flattening in cores with stars in our samples. However, it is possible 
that central flattening was not observed due to the limited spatial resolution (~ 18 arcsec or 
2500 AU) and data sampling (~ 21 arcsec or 2900 AU) in our observations. Our result shows 
that cores with stars have shallow power-law density (or integrated intensity) distribution. 
It is interesting that this is in contrast to the result in the lower spatial resolution study by 
Caselli et al. (2002c), where cores with stars show steep integrated intensity distribution. 
There is a possibility that the shallow density profile in cores with stars can be results of the 
core collapse (Foster & Chevalier 1993) and/or core dispersal due to the molecular outflow. 
The present observations are probably not enough for disentangling these possibilities, and 
further observations (high-resolution dust continuum map, near-infrared color-excess map, 
higher resolution N 2 H + imaging, etc) are desirable. 



5. Summary 

On the basis of N 2 H + observations toward Taurus, we have studied the physical prop- 
erties of the molecular cloud core. The core radius, linewidth, and intensity distribution 
do not much differ between starless cores and cores with stars. This result is in contrast 
with that previously obtained in H 13 CO + . We suggest that depletion of H 13 CO + causes this 
difference. From the critical pressure analysis for Taurus cores, there is no systematic differ- 
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ence between starless cores and cores with stars. Both are not far from the critical state for 
equilibrium. We suggest that the starless cores which are almost thermally supported evolve 
toward star formation by keeping close to the critical state. This result is in contrast with 
that obtained in the intermediate-mass star forming region OMC-2/3, where the molecular 
cloud core evolves by dissipating turbulence largely. The density profile is investigated from 
the integrated intensity distribution in the cores. Cores with stars show shallow density 
profiles, r~ L8 to r _L6 . 

K. T. is grateful to Takenori Nakano for comments on the draft and to Jeong-Eun Lee 
for discussion. 
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Fig. 1. — Two-dimensional distribution of the velocity- integrated intensity of the main N 2 H + 
J = 1 — > component (-F\, F = 2, 3— >1, 2) toward eight molecular cloud cores. The lowest 
contour level and level interval are 0.05 K km s -1 except for core L1551, for which they are 
0.1 K km s -1 . Small dots represent the observed positions. Crosses represent the positions 
of young stars. 




Fig. 2. — Composite N 2 H + spectra made by using the central nine positions. The result of 
the hyperfine fitting is also shown. 
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Fig. 3. — The hyperfine fitted linewidth of N2H 4 " is plotted against the angular distance from 
the core center (impact parameter) b. The open circle and error bar represent the average 
and la deviation of the linewidth at the same angular distance, respectively. 



- 15 - 



10 



8 



o OMC-2/3 w/o STAR 
• OMC-2/3 w/ STAR 
□ Tau w/o STAR 
■ Tau w/ STAR 
S 



'E 
o 



o 
CL 



10 { 



Ps(Tau) 



10* 



10' 



0.1 



SI 



□ 



I 



n 1 — i — i — i r~_ 



o 



Q 



Ps(OMC-2/3)- 



o 



_i i •. i i i 1 1 



_1 I I I ! I l_ 



1 10 

CORE MASS M (Mo) 



100 



Fig. 4. — The critical pressure P cr is plotted against the core mass M. For Taurus cores with 
stars, P cr corrected for mass lost due to accretion onto the protostar and due to the outflow 
is also shown as the filled triangle. For reference, OMC-2/3 cores (Aso et al. 2000) are also 
plotted. The horizontal lines at 4 x 10 5 K cm -3 and at 2 x 10 6 K cm~ 3 represent the estimated 
core surface (external) pressure P s for Taurus cores and OMC-2/3 cores, respectively. 




Fig. 5. — The velocity-integrated intensity of the main N 2 H + J = 1 — > component (Fi, F 
= 2, 3— >1, 2) plotted against the angular distance from the core center (impact parameter) 
b. The filled circle and error bar represent the average and la deviation at the same angular 
distance, respectively. The intensity profile is fitted with a power law I oc b~ p convolved with 
the telescope beam. For Mizl and Miz2, the fitting results by excluding two inner radius 
bins are also shown as broken lines. 
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Table 1. Observed molecular cloud cores. 



Field 


R.A.(J2000) a 


Dccl.(J2000) a 


Parent cloud name 


IRAS/young star 


Class 


Mizl 


4 18 10.3 


27 36 03 


L1495 






Miz2 


4 19 37.0 


27 15 35 


L1495 






L1521F 


4 28 39.3 


26 51 43 


L1521 






Miz7 


4 19 41.5 


27 13 32 


L1495 


04166+2706 


I 


Miz8 


4 19 58.2 


27 10 22 


L1495 


04169+2702 


I 


Miz8b 


4 19 52.2 


27 11 45 


L1495 






L1551 


4 31 32.5 


18 08 25 


L1551 


L1551 IRS5 


I 










L1551 NE 





L1527 


4 39 53.2 


26 03 47 


L1527 


04368+2557 






a The maximum position of the main N2H+ J = 1 — > component F = 2, 3— >1, 2) is 

listed. Units of right ascension are hours, minutes, and seconds, and units of declinations are 
degrees, arcminutes, and arcseconds. 
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Table 2. Radius, velocity, linewidth, optical depth, and excitation temperature. 



Core 


R (pc) 


V LSR (km s-!) b 


Ad (km s" 1 ) 


TOT b 


T ex (K) b 


Mizl 


0.038 


6.75±0.01 


0.28±0.04 


d 


d 






(6.76±0.01) 




(2.81±1.07) 


(4.86±0.62) 


Miz2 


0.036 


6.83±0.01 


0.24±0.05 


4.00±1.33 


6.23±0.83 






(6.84±0.01) 




(2.73±0.86) 


(5.99±0.80) 


L1521F 


0.029 


6.50±0.01 


0.28±0.04 


9.72±1.79 


5.07±0.23 






(6.51±0.01) 




(9.63±0.87) 


(4.65±0.09) 


Miz7 


0.024 


6.72±0.01 


0.24±0.05 


7.64±1.33 


5.78±0.31 






(6.70±0.01) 




(5.03±0.72) 


(5.47±0.26 


Miz8 


0.029 


6.75±0.01 


0.33±0.05 


7.92±1.78 


5.02±0.29 






(6.77±0.01) 




(4.54±1.06) 


(4.78±0.33) 


Miz8b 


0.036 


6.73±0.01 


0.23±0.06 


3.96±1.79 


5.91±1.03 






(6.71±0.01) 




(1.89±0.82) 


(7.40±1.67) 


L1551 


0.038 


6.41±0.02 


0.40±0.10 


1.06±0.79 


11.20±5.53 






(6.43±0.01) 




(2.19±0.42) 


(7.60±0.73) 


L1527 


0.032 


5.85±0.01 


0.27±0.07 


7.43±1.94 


4.20±0.22 






(5.88±0.01) 




(4.21±0.88) 


(4.60±0.28) 



b Upper row shows the value at the intensity peak, and the lower row in parentheses 
is the value from the composite spectrum of central nine positions. The error represents 
lcr error in the hyperfine fitting. 

intrinsic N2H+ linewidth. For Miz8, L1551, and L1527, the values derived for the 
core edge (b > 40", 80", and 40", respectively) are used to minimize the effect of the 
protostellar collapse and molecular outflow. For the remaining sources, the average values 
over the cores are used. 

d Very large uncertainty. 
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Table 3. Column density, mass, critical pressure, radial velocity distribution, and radial 

intensity distribution. 



Core 


AT(H 2 ) (cm- 3 ) c 


M (M ) f 


Per (10 5 K cm- 3 )e 


Power index of linewidth h 


Power index of intensity distribution 1 


Mizl 


(1.0 ± 0.4 x 10 22 ) 


1.1 


7.3 


-0.23±0.18 


(C.C.=0.88) 


-1.2±0.3 k 


Miz2 


1.3 ± 0.5 x 10 22 
(1.1 ± 0.4 x 10 22 ) 


1.0 


6.1 


-0.42±0.14 


(C.C.=0.79) 


-1.1 ± 0.2 k 


L1521F 


3.0 ±0.6 x 10 22 
(2.8 ± 0.3 x 10 22 ) 


2.2 


1.6 


-0.17±0.08 


(C.C.=0.48) 


-1.0 ±0.2 


Miz7 


2.2 ± 0.4 x 10 22 
(1.4 ± 0.2 x 10 22 ) 


1.0 


3.1-6.4 


-0.02±0.34 


(C.C.=0.98) 


-0.7 ± 0.2 


Miz8 


2.6 ± 0.6 x 10 22 
(1.7 ± 0.4 x 10 22 ) 


1.3 


4.4-7.5 


-0.41±0.36 


(C.C.=0.92) 


-0.6 ±0.1 


Miz8b 


1.3 ± 0.6 x 10 22 
(0.9 ± 0.4 x 10 22 ) 


0.9 


6.7 


-0.14±0.25 


(C.C.=0.68) 


-0.8 ±0.1 


L1551 


3.2 ± 2.7 x 10 22 
(3.0 ± 0.6 x 10 22 ) 


3.0 


1.5-2.6 


-0.31±0.13 


(C.C.=0.97) 


-0.8 ±0.1 


L1527 


2.2 ± 0.6 x 10 22 
(1.4 ± 0.3 x 10 22 ) 


1.1 


5.0-7.1 


-0.23±0.19 


(C.C.=0.79) 


-0.7± 0.1 



°Upper row shows the value at the intensity peak, and the lower row in parentheses is the value from the composite spectrum of central 
nine positions. The error represents 1<t error in the hyperfine fitting. 

'Mass lost due to the forming star and the associated molecular outflow is not taken into account. 

s For cores with stars, the "range" represents the values uncorrected (upper bound) and corrected (lower bound) for the mass loss due 
to accretion onto the star and due to the molecular outflow. 

h Power-law index q for the linewidth-size relation (Fig. 3) estimated by using the non-linear least-squares fitting in the form Av oc b q . 
The correlation coefficient is listed in parentheses. 

'Power-law index from the fit to the radial integrated intensity distribution of the main component (Fig. 5). 

JVery large uncertainty. 

k Central part is excluded for fitting. 



